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ABSTRACT 


A thermal/fluids analysis of a direct gain solar thermal upper stage engine is presented and the 
results are discussed. The engine was designed and constructed at the NASA Marshall Space 
Flight Center for ground testing in a facility that can provide about 10 kilowatts of concentrated 
solar energy to the engine. The engine transfers energy to a coolant (hydrogen) that is heated and 
accelerated through a nozzle to produce thrust. For the nominal design values and a hydrogen 
flowrate of 2 lb./hr., the results of the analysis show that the hydrogen temperature in the 
chamber (nozzle entrance) reaches about 3800°F after 30 minutes of heating and about 3850°F at 
steady-state (slightly below the desired design temperature of about 4100°F). Sensitivity 
analyses showed these results to be relatively insensitive to the values used for the absorber 
surface infrared emissivity and the convection coefficient within the cooling ducts but very 
sensitive to the hydrogen flowrate. Decreasing the hydrogen flowrate to 1 lb./hr. increases the 
hydrogen steady-state chamber temperature to about 4700°F, but also of course causes a 
decrease in thrust. 
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INTRODUCTION 


Other means of producing thrust than by chemical propulsion are being investigated in efforts 
to find better and cheaper ways to explore space and insert payloads into Earth orbit. One of 
these innovative approaches is to use solar energy to produce thrust. A solar thermal engine 
collects concentrated solar energy which is used to heat a propellant to temperatures on the order 
of 4000°F to 5000°F. The heated propellant exits at a high velocity through the engine’s nozzle, 
thus producing thrust. There are technology limitations to these concepts though. Absorber 
materials such as tungsten, rhenium or molybdenum capable of withstanding these high 
temperatures are required. Also, flexible concentrators large enough to focus enough energy into 
the absorber with accurate pointing and stiffness along with being lightweight are a considerable 
technical challenge. 

The major benefit to using the solar thermal engine is the expected increased specific impulse 
(I sp ). This high I sp is due to the low molecular weight of H 2 and the exhaust velocities of the 
propellant exiting the nozzle. The I sp is proportional to \JTJ M , where T c is the propellant 
stagnation temperature and M is the propellant molecular weight [1]. Since H 2 has the lowest 
molecular weight of any element, increasing the temperature is necessary to increase the I sp . 
Typical L0 2 /LH 2 chemical stages have I sp of 436 sec., whereas a typical Solar Thermal Upper 
Stage (STUS) design should have an I sp of about 860 sec. if a fluid temperature of approximately 
4100°F could be attained. The STUS, though, is designed to produce about 2 lb. of thrust, which 
only makes it useful to boost payloads to higher orbits over a long period of time. 
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Using solar energy in a propulsion system was first proposed by Kraft Eriche in 1956[2], No 
other work was done in this field until the early 1960’s. During this time the Air Force’s Rocket 
Propulsion Lab through a contract with Electro-Optical Systems, Inc. performed detailed 
analyses of the solar thermal rocket concept and carried out subscaled demonstrations of solar 
heating of the propellant. A Solar Hydrogen Rocket was developed. A five-foot reflector 
focused solar energy onto a coil through which hydrogen flowed. The solar absorbers were made 
of molybdenum and rhenium. Hydrogen temperatures were measured up to about 4,400°F which 
corresponded to a calculated I sp of about 700 sec. 

Several solar absorber/thruster concepts have been examined. Some of the designs have 
windows which allow the solar energy into the engine. In these designs refractory seedant 
material is used to absorb the energy and transfer heat to the working fluid. The seedant material 
and the fluid would be expelled from the engine’s nozzle. A windowless concept has the solar 
energy being absorbed by the interior walls of the absorber/thruster and the heat is collected when 
the working fluid flows over the surface of the absorber. Though the windowed designs may be 
capable of higher temperatures, the windowless (direct gain) design is much simpler and was 
found to meet performance requirements [2], 

In 1995 work began at NASA’s Marshall Space Flight Center (MSFC) on a Center Director’s 
Discretionary Fund (CDDF) Project to build and test a direct gain solar thermal engine. This 
engine would use hydrogen (H 2 ) as its propellant. The work presented here focuses on the 
thermal/fluids modeling of the engine for the NASA CDDF Project. The models developed were 
used to calculate the pressure loss in the engine, fluid velocities, and the fluid temperature along 
the solar thermal absorber up to the nozzle. Using this information one can determine the 
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performance of the engine. A sensitivity study was performed to investigate the effects of 
changes to several of the parameters that influence the engine’s performance: propellant flow 
rate, convection coefficient, and absorber optical properties. 

THE SOLAR THERMAL ENGINE 


Design 

The solar thermal engine is shown in Figure 1 . The absorber cavity was made of tungsten 
which has a melting point of 5,684 °F. The absorber wall was approximately 0.08 in. thick. The 
absorber cavity, which totaled 16.3 inches in length, was a cylinder with a hemispherical end. 
The diameter at the mouth of the absorber was 2.652 in. while the diameter at the end of the 
cylinder was 2.588 in. The cylindrical portion of the absorber was 15 in. long with the balance 
being the end cap. The tapering of the cylinder was to aid in the removal of the absorber from the 
mandrel on which it was sprayed using Vacuum Plasma Spray (VPS) techniques. The H 2 ducts 
were formed using VPS techniques also. Tungsten was deposited on a mandrel that was threaded 
such that it produced a helical duct when the mandrel was removed. There were two parallel 
ducts each having a helical pitch of 1 , meaning that one helical revolution of the absorber was 
completed for every inch of the absorber length. The shell that formed the ducts had an outer 
diameter of 3.213 in. at the beginning and a diameter of 3.170 in. at the end. The ducts were 
triangular in cross section with an area of 0.045 in 2 ., base of 0.45 in. and height of 0.21 in. The 
walls of the ducts were 0.08 in. thick. There was a gap of 0.204 in. between tip of the duct walls 
and the insulation. The insulation was Fiberform®, a high temperature insulation that was 
designed to operate at temperatures up to 5000°F. The Fiberform® insulation facing the ducts 
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had a low emittance graphite coating, providing a reflective coating effective at high temperature 
operation. The collar at the mouth of the absorber was made of nickel. The entire engine and 
insulation assembly was enclosed in stainless steel. 

The energy used to heat the engine is provided through concentrated solar energy. The solar 
energy spectral distribution is approximately 7 percent ultraviolet, 46% visible and 47% near 
infrared. The majority of the energy emitted by the sun is at a much shorter wavelength than the 
infrared energy emitted by a body near the design temperature of the surface of the absorber. 
This allows for the selection of finishes that are absorbing or reflective in the solar spectrum 
while performing oppositely at long infrared wavelengths. For an engine of this type, it is 
desirable to have a high solar absorbtivity and a low emissivity at the temperature of the surface, 
because that would enable more energy to be transferred to the fluid rather than be radiated to the 
environment. The absorbitivity (a) is the fraction of energy that is absorbed by a surface 
compared to the energy absorbed by black surface (a=l), while the emissivity (e) is the fraction 
that is emitted due to the temperature of the surface compared to a black (e=l) surface of the 
same temperature. The design specification of the engine required that the surface of the first 4 
inches of the absorber cavity be polished to have an a of 0.25 while the balance of the absorber 
would have an a of 0.75 at solar radiation wavelengths. The a of the first 4 inches of the 
absorber was reduced to help increase the reflectance in that area. 

Thermal Model 


6 


A thermal model was developed to represent the conditions in which the engine would be 
tested. The engine would be tested in a vacuum chamber; thus it was assumed that the heat 
transfer from the outer surface of the engine would be due only to radiation. Also the solar flux 
input is based on what is expected to be produced in the test environment. No attempt was 
made to model the nozzle flow or heat transfer in the nozzle area. Only the heat transfer and 
fluid dynamics associated with the absorber, ducts and the insulation were modeled. A detailed 
description of the model and the results of the analysis are given in reference 3. 

The solar thermal engine was modeled using SINDA (Systems Improved Numerical 
Differencing Analyzer) [4], The SINDA model was a set of nodes that were used to discretize 
the different components of the engine. A network of nodes simulating all parts of the engine 
was connected using conductors: radiation, convection, and conduction. SINDA uses finite 
difference techniques to solve the network for nodal temperatures. 

Radiation view factors were calculated using TRASYS (Thermal Radiation Analysis System) 
[5], which is a program that has the general capability to solve the radiation aspects of a thermal 
analysis problem. It can be used to calculate the intemode radiation interchange of a model as 
well as calculate the incident and absorbed heat rate data from environmental heat sources. For 
this analysis it was used to solve the complex radiation interchange between all of the nodes in 
the absorber and the environment. TRASYS was also used to calculate the radiation exchange 
between the flow channel walls and the insulation. 

Figure 2 is a schematic of the absorber nodes as they would appear if the cylinder has been 
sliced down one side and unrolled. The absorber cavity was divided into 15 rings, which were 1 
in. in width. The rings were then each divided circumferentially into 12 equal segments. This 
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formed 180 nodes that made up the cylinder of the absorber cavity. The endcap was a single 
node that was a 1.3 in radius hemisphere. Figure 3 shows the H 2 flow paths over the absorber 
elements. There were two separate parallel flow paths (ducts), each of width 0.45 in. at the base 
of the duct cross-section. 

Fluid Flow Model 

Of particular interest in modeling the engine was the temperature of the fluid as it flowed 
along and exited the absorber. The primary mode of heat transfer to the fluid was forced 
convection, so part of the fluids modeling focused on determining the convection coefficients 
needed to calculate the temperature of the fluid. One hundred and eighty fluid nodes per duct 
were used to simulate the H 2 flowing through the engine. Single-phase (gas) fluid flow was 
assumed since this is the fluid condition that is expected during testing of the engine. The 
pressure loss of the fluid as it flowed through the engine ducts was also of interest. Both of these 
were complicated by the shape of the duct and also curvature effects due to the coiled 
configuration of the ducts. As described in detail in reference 3, these effects were taken into 
account in the model, and the curvature had a significant effect on the value of the convection 
coefficient in the ducts. 

Incident Solar Energy Flux Distribution 

The total spectrum has to be considered when examining the radiation environment for the 
solar thermal engine: the visible for solar energy input and the infrared for radiation of absorbed 
energy. Since the heat input to the absorber is solar energy (visible light) the optical properties 
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for the calculations were values that were integrated over the solar spectrum, for the tungsten 
surface. QFLUX, a Monte Carlo simulation tool, was used to determine the light distribution 
inside the absorber cavity used for the solar energy input for the analysis. The code divided the 
cavity into small nodes then generated random light rays entering the cavity. The light flux 
distribution inside the cavity was calculated by determining where each of the rays struck the 
surface of the absorber and where its reflections struck. The code calculated how much of the 
energy was absorbed based on the solar absorbtivity of the surface, and the balance of the energy 
was reflected. The light ray was continually reflected until its intensity was only 1% of the 
original value. 

The output from the simulation was a normalized distribution of the light on the absorber 
specific to the optical properties of the absorber and the angle at which it was concentrated. In 
this case the light had been concentrated such that it had a half angle of 32° as it entered the 
absorber cavity. The flux data was given in terms of R, where R is the amount of energy 
deposited on a single absorber element in the simulation and is defined by 


R = 




(i) 


where q xy is the light flux over the area a xy and Q T is the total energy passing through the area of 
the mouth A mouth of the absorber. 

The light distribution profile is shown in Figure 4. The nominal design was specified to 
have solar wavelength absorbtivity of 0.25 for the initial 4 (polished) inches of the absorber and 
0.75 for the remainder, with corresponding infrared wavelength emissivities of 0.1 and 0.22. The 
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distribution has two peaks at about 2.5 in. and 5 in. along the absorber cylinder. The first peak is 
due to the light entering the absorber at a half cone angle of 32° and striking the surface and being 
absorbed. The second peak is caused by the reflections of light from the first region of incidence. 
Light is reflected further into the absorber, but it is more evenly distributed on the absorber 
surface. The peak at the end of the absorber is due to light directly incident on the surface from 
the source. 


RESULTS 

The cases analyzed were based on tests planned to be performed at the Air Force Research 
Laboratory (AFRL) test facility. Ten kilowatts of energy was specified to be concentrated into 
the absorber cavity of the engine, with gaseous H 2 flowing through the absorber at the rate of 2 
lb./hr for the nominal design. The engine was assumed to be in a vacuum chamber and at an initial 
temperature of 70°F. (For all cases analyzed, the results showed that the flow was laminar, the 
Mach number was less than 0.1, and the pressure drop from the duct inlet to the chamber was 5 
psi or less.) 

Shown in Figure 5 are steady state results for the fluid temperature distribution along the 
length of the absorber, with the fluid chamber temperature reaching about 3850 °F. The transient 
fluid temperature distribution is shown in Figure 6, where it is seen that after 30 minutes of 
operation the fluid chamber temperature is within about 40 °F of the steady state temperature 
and that even after only 5 minutes of operation the temperature is within 500 °F of steady state. 
Steady state axial distributions of temperatures at different radii at the 0° circumferential location 
are presented in Figure 7, and the dominating role of the heat transfer resistance of the insulation 
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layer is quite evident. Results of a study in which the insulation thickness of the nominal design 
was varied are shown for axial temperature distributions in Figure 8 and for the fluid chamber 
temperature in Figure 9. It appears relatively little is gained in terms of fluid chamber 
temperature by increasing the insulation thickness above about 0.6 inch. 

Additional cases were run to investigate the sensitivity of the results from the model to 
several of the input variables. It was found that increasing or decreasing the convection 
coefficient in the ducts by 50% had a negligible effect on the temperature of the fluid at the 
chamber. 

Cases were also run in which the emissivity of the absorber surface was increased or 
decreased by 50%. The results of the radiation heat transfer sensitivity analysis show that with 
lower emissivity values the temperature of the H 2 increases. Less heat was radiated to the 
boundary through the mouth of the absorber, which means that the fluid received more heat and 
experienced increased temperatures. 

Table 1 shows that when emissivity was reduced by 50% the amount of heat that was 
radiated through the mouth of the absorber was 2640 BTU/hr., and the fluid exit temperature was 
3904°F. In the nominal case the radiated heat was 3104 BTU/hr. and the fluid temperature was 
3847 °F. And finally, in the case where emissivity was increased by 50% the radiated heat was 
341 1 BTU/hr., and the exit temperature was 3809 °F. Table 1 shows the relative amount of heat 
that was radiated by the absorber as the emissivity of the absorber was changed. Nine percent of 
the total heat (lOkW) was radiated for the nominal case whereas only 7.7% was radiated for the 
50% reduced case, while 10% was radiated for the 50% increased case. This table also shows the 
amount of heat that was radiated to the environment from the outer shell of the engine. The 
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engine performance is thus relatively insensitive to e. Multiplying e by a factor of 3 only 
changed the chamber temperature from 3900°F to 3800°F. 

The fluid chamber temperature is very sensitive to the hydrogen flow rate, as shown in Figure 
10. Decreasing the flow rate from 2 Ib./hr. to 1 lb./hr increased the chamber temperature to about 
4700 °F. Of course, the reduced mass flow rate also would reduce the thrust. 

SUMMARY AND CONCLUSIONS 

A direct gain solar thermal engine has been analyzed and the results presented. Shown was 
that the STUS engine design, with input energy of lOkW, will have a chamber temperature of 
3850°F which corresponds to an I sp of about 820 seconds. The Ap through the engine was 
calculated to be about 4.9 psi. 

Sensitivity analyses show that the major factor that affects the temperature of the fluid is the 
mass flow rate of the working fluid if the collected solar energy is held constant. Reducing the 
hydrogen mass flow rate to 1 lb./hr. increases the chamber temperature to about 4700°F, but with 
a corresponding decrease in thrust. 

The temperature of the fluid is sensitive to the convection coefficient (h), but not to a large 
extent. For h to have a major effect on the temperature, the convection coefficient would have to 
be changed by a factor greater than 2. This implies that a simple geometry change of the duct is 
not enough to increase the fluid temperature. Similarly the emissivity of the absorber surface 
must be much lower to increase the chamber temperature of the fluid. Also implied in these 
sensitivity analysis results is that the uncertainty in the emissivity and the fluid properties used 
to calculate h has a minimal effect on the results. 
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When developing a flight engine the mass and volume occupied are also very important. 
Reducing the insulation thickness of the engine reduces the mass and volume but also reduces the 
engine’s chamber temperature while increasing the outside shell temperature. Results of the 
insulation sensitivity analysis show that with insulation approximately 0.6 in. thick the steady- 
state chamber temperature will be approximately 3740°F while the stainless steel outside shell 
temperature will be about 1000°F. This chamber temperature is only about 100°F less than the 
chamber temperature of the engine with 3.2 in. of insulation, which has four times the volume of 
an engine with 0.6 in. of insulation. 
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Figure 4. Incident Light Flux Distribution 



Figure 5. Fluid Temperature along Absorber Length 
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Figure 6. Axial Temperature Distributions - Transient Results 
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Figure 7. Axial Distributions of Temperature at Different Radii 





Figure 8. Axial Temperature Distribution of Outer Shell as a Function of Insulation 
Thickness 
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Figure 9. Fluid Chamber Temperature as a Function of Insulation Thickness 
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Figure 10. Influence of Hydrogen Flowrate on Fluid Temperature 



Table 1. Heat Transfer Results from Radiation Sensitivity Analysis 


Factor Applied to 
£ 

Heat Reradiated 

Through 

Absorber Mouth 

(BTU/hr.) 

Heat Radiated 
from Outer Shell 
(BTU/hr.) 

Total (BTU/hr.) 

Chamber 
Temperature (°F) 

0.5 

2641 

1218 

3859 

3904 

1.0 

3104 

1311 

4414 

3847 

1.5 

3411 

1373 

4784 

3806 


Table 2. Absorber Heat Balance 































